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Abstract 

We present the Yano-Koonin-Podgoretskii parametrisation of the correlation 
function. Compared to the conventionally used Cartesian parametrisation, this 
one provides more straightforward measurement of the duration of the emission 
process in the fireball and a clearer signal of the longitudinal expansion, which is 
expected in ulrarelativistic heavy ion collisions. 
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1 Introduction 



The study of ultrarelativistic heavy ion collisions is motivated by the prediction of QCD 
lattice gauge theory that hadronic matter undergoes at energy density of 1 GeV/fm a 
phase transition into the quark-gluon plasma (QGP). To determine the energy density 
obtained in heavy ion collisions, a measurement of the dimension and lifetime of the 
fireball is needed. However, due to very short lifetime of the examined object it is 
impossible to use conventional methods as e.g. the scattering of external particles. 
The most direct measurement of spatio-temporal characteristics of the collision region 
is provided by Hanbury-Brown/Twiss interferometry, a method developed originally in 
radioastronomy in the fifties pj|. 

In the field of particle physics the idea was first applied by Goldhaber, Goldhaber, 
Lee and Pais in 1960 g. 

In the last years this method has been widely used in ultrarelativistic heavy ion 
collisions to investigate the spatio-temporal characteristics of the fireball. Important 
conceptual extensions have been connected with this application. This is mapped by 
numerous reviews and introductory articles in the field ||, ^j, pi ^, m] . In this article we 
review a new parametrisation of the correlation function, and point to some of its most 
important features. We assume complete chaoticity of the source. This assumption has 
not been proven yet. However, it might be supported by the successful experimental 
use of this method. 

In the following Section we start by recalling a few basic notions and relationships 
of the HBT interferometry. We pay special attention to the problems arising from the 
on-shell constraint for the bosons in the final states and from the shape of the correla- 
tion function, which is very similar to Gaussian even for many different particle sources. 
We recall the way how to handle these problems. This leads to the parametrisation of 
the correlation function by a Gaussian function in three dimensions. Sections || and 
H] are dedicated to these Gaussian parametrisations and to the spatio-temporal inter- 
pretation of their parameters. As a first step in Section || we shortly recall the basic 
properties of the conventional, so-called Cartesian parametrisation. Then we present 
in Section || the recently derived Yano-Koonin-Podgoretskh (YKP) parametrisation. 
We introduce its properties, pay particular attention to the interpretation of the pa- 
rameters and we compare the quality of the information about the source obtained via 
this parametrisation with the information from the Cartesian one. The advantageous 
features of the YKP parametrisation reside especially in the direct measurement of 
the emission duration in the fireball and a straightforward signal of the longitudinal 
expansion. We also comment shortly on the technical difficulties which can arise in 
fitting the data to this parametrisation. 



2 Correlation basics 



We start the brief review of basic theory with the formula describing the connection 
between the source and the observed correlations (3|, ||, ^, 10 1 



C(q,K) 



\Jd 4 xe ic >- x S(x,K)\ 
\JSxS(x,K)\ 2 



(1) 
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Here C(q, K) is the correlation function depending on the variables of average momen- 
tum K = \{k\ + fo) and the momentum difference q = k\ — fo- The emission function 
S(x,p) expresses the particle source. It is the phase-space (Wigner) density describing 
the probability for the boson with four-momentum p to be produced of the space-time 



point x [§, p|, |10|. 



In the experiment the correlation function is measured, and one would like to obtain 
the information about the source, i.e., to extract S(x, K) from the correlation function 
C(q,K). Unfortunately, the measured information about S(xK,) is not unambiguous 
due to the on-shell constraint for the particles in the final state. The momenta fulfill 
the relation 

q ■ K = (2) 

which can also be written in the form 

q ° = P- q , A = |L«|±. (3) 

Due to this relation, eq.(Q) is not invertible, because only three of the four components 
of the momentum difference are independent. Hence, the analysis of correlation data 
is not possible in a completely model independent way, and model studies have to be 
employed. 

Another point which we want to recall here is connected to the fact that the cor- 
relation function has in the region of small q a Gaussian shape for a very wide class of 
source functions [Q, J5|, at least for thermal models^. Then the emission function can 



be written in the form |6|, |17J , [L8|, [lg [20| 

' +5S(x,K). (4) 



-^(K)iV(K)r(K) 



S(x, K) = N(K) S(x(K),K) exp 

Here the space-time coordinates relative to the "source center" x(K) 

x » = x » - #*(K), x^(K) = (x 1 *) (5) 
with the space-time averages over the source defined as 

The statement is that only the second moments given by the matrix B^ v are measured. 
The first two terms in the Gaussian part on the r.h.s. of eq.(Q) are not seen in the 
correlation measurement due to the normalization to the single-particle spectra in 
eq.(||). The correction to the Gaussian approximation 5S(x,K) is neglected as far as 
one is fitting the data by a Gaussian. For the class of models to be discussed below, 



non-gaussian effects from 5S have been shown to be small [17]. 
Inserting (^|) into (|l]) and neglecting 5S one gets 



C(q,K) = l + e ~rr^ *V . (7) 



1 Recently it has been found that the Gaussian shape of the correlation function is not preserved if 
some of the bosons originate from resonance decays |n], ^| ^] . 
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This shows that the inverse of the matrix B^ u , expressible as 



(8) 



are measured. However, due to the on-shell constraint, not all 10 components of this 
matrix are measured. It is possible to measure only 6 combinations of the matrix 
elements. In what follows, we will treat the case of azimuthally symmetric sources. 
This corresponds to central collisions. In this case 7 non- vanishing matrix elements 
are needed for the full description of the source according to eq. @j , but the number of 
measurable combinations of them is reduced to 4. 

The question investigated in this work is how to choose the Gaussian parametrisa- 
tion of the correlator such that the interpretation of the four measured parameters in 
terms of the space-time dimensions of the source is as simple (and model-independent) 
as possible. 

Interpreting the correlation parameters (called also HBT radii) it is important 
to realize that the correlations do not measure the dimensions of the whole source. 
According to eqs.(0) and (||) they are combinations of the space-time variances of the 



source and they are measuring only the so-called homogeneity lengths [16, 17, 21, 22 



23]. It is intuitively obvious that if the measurement is focussed on bosons with given 
average momentum only the part of the source producing these bosons is measured. 
This is relevant e.g. in sources with strong expansion. Then a particular momentum 
is not produced by the whole fireball, but only by a small part of the collision region. 
Bosons with different momenta originate from different parts of the source. Hence the 
HBT radii can be momentum dependent, and this dependence can carry information 
about the dynamics of the source. It is therefore an object of very intensive theoretical 



[17, 21, 24, 25] and experimental [26, El 28, 29 1 study. 



3 Gaussian parametrisations 1: 
Cartesian parametrisation 

In this Section we shortly recall the basic properties of the Cartesian parametrisa- 
tion. This sets the stage for a comparison to the Yano-Koonin-Podgoretskii (YKP) 
parametrisation in the next Section. 

We use the usual coordinate system with the z-axis in the direction of the beam and 
the x-axis parallel to the transverse component of the average momentum K. Then 
the z-axis is called also longitudinal, the x-axis is labeled as outward and the remaining 
y-direction is denoted as side-ward. 

Both Cartesian and YKP parametrisations are constructed for azimuthally sym- 
metric events. As already mentioned, in this case four fit parameters are needed for 
the Gaussian parametrisation of the correlation function. Different parametrisations 
can be obtained from eq.(|7|) by different choices for the three independent components 
of the momentum difference q according to eq.(||). 



In the "standard" Cartesian parametrisation |16|, 22] the three independent q— 
components are chosen to be the three spatial components of the momentum difference. 
The parametrisation is given by the following formula 

C(q,K) - 1 = exp[-i^(K)g s 2 - R 2 (K)q 2 - R 2 (K)q 2 - 2 R 2 ol (K)q sqi ]. (9) 
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The interpretation of the four HBT radii appearing in this parametrisation is most 
straightforward in the so-called LCMS (Longitudinal Co-Moving System) frame. This 
is the longitudinally boosted frame in which the average pair momentum has only a 
transverse component, K[ = 0. The interpretation of the HBT radii in this frame is 



follows from the model independent expressions [16, 3C] 



R 2 — 


<« 2 > 


(10) 


R 2 Q = 


{(x - 0jf) 


(11) 


Rf = 


(?) 


(12) 


R 2 oi = 


((x-(3j)z) 


(13) 



It is clearly seen that these radii mix spatial and temporal information about the 
source. It has been found, that the emission duration can be measured by the difference 

Rl - R 2 S = Pi (?) - 2/3jl (xi) + (x 2 - y 2 ) . (14) 

In the typical case of ultrarelativistic heavy ion collisions the last two terms at the r.h.s. 
of this equation can be treated as perturbation [JDJ and the difference is a measure for 
the effective emission duration (the so-called lifetime of the source). However, two 
problems are connected with this measurement. The first is caused by the presence of 
the pre- factor (3 2 _, which makes this observable small in the region of most data points. 
The second one is more serious. In typical measurement in ultrarelativistic heavy ion 
collisions R 2 and R 2 are typically bigger than the expected lifetime. Thus the lifetime 
measured in this way is obtained with big statistical errors. 



4 Gaussian parametrisations 2: 

the Yano-Koonin-Podgoretskii parametrisation 

The Yano-Koonin-Podgoretskii parametrisation |32| , [33] , |34|| is given by the following 
formula 

<7(q,K)-l = exp[-tfi(K) g i-i?^ (15) 

In this parametrisation, the three independent momentum difference components are: 
<7_L = \A?o + Q 2 i Qi an d Q°- The fit parameters are now R±, R 2 , Rq and the so-called 
Yano-Koonin velocity v appearing here in the four-velocity U. This four velocity is 
assumed to have only a longitudinal component 

U(K) = 7 (K) (1, 0, ,0, v(K)), 7 (K) = 1 . (16) 

y/1 - V 2 {K) 

The YKP parametrisation is constructed in the way that the three radius parame- 
ters are invariant under longitudinal boosts, i.e., in any longitudinally boosted reference 
frame the analysis of the data should give the same results. 



4 



The model independent expressions for the HBT radii have been found in refs.[18 
[l9tl • They are easiest written with the help of notational shorthands 



A 



B 



C 




(17) 
(18) 
(19) 



Here £ = x + iy. Furthermore we use that for azimuthally symmetric sources (y) = 
(xy) = and (£ 2 ) = (x 2 — y 2 ). Then the HBT radii and the YK velocity are given by 
the following formulae 



Rl 
Rj 

Ri 




2C 
A + B 



(20) 
(21) 

(22) 
(23) 



Since eqs.(^) and (15) are just different parametrisations of the same correlation 
function, there must exist a simple connection between the YKP parameters and those 
of the Cartesian parametrisation. It is given by [18] 



r?2 — r?2 



R 2 . 



R 



R-li 



R ± , 

Ph 2 (Rl + v 2 R\ 
/3f)4 + 7 2 (A 



0± {-m\ 



,f (rI + r\ 



7 2 (A 



v) [Rq + R\ 



(24) 
(25) 

(26) 

(27) 



Although the three radius parameters are longitudinally boost-invariant, there is 
a special frame in which their interpretation is simplest. The frame is defined by the 
relation v = and it is called the Yano-Koonin (YK) frame. In this frame the formulae 
dtj-ED simplify to 



Rl = 


(0 


5 


R 2 = 


<*> 




Ro = 


(?) 


2 / ~ ~\ 1 / 



X 



(28) 
(29) 

(30) 



where x, y, z, t are now measured in the YK frame. It has been pointed out in [18] 
and shown in a detailed numerical model study in [19] that the last two terms on the 



r.h.s. of eqs.(29) and (|30| ) can be considered as a small correction within a class of 
thermal models with Gaussian geometrical density profile. Hence, R» measures (up 
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to the corrections) the dimension in the longitudinal direction in the YK frame and 
i?o is a measure for the emission duration in that frame. In contrast to the Cartesian 
parametrisation, the time variance appears here as the leading term of the fit parameter 
Rq. Thus the problems with the accuracy arising by the subtraction of two numbers 
of the same order are avoided here. 

The YK velocity v has been found to coincide up to small corrections with the 
longitudinal velocity of the point of maximal emissivity for bosons of given momentum. 
This will be illustrated in the model study in Subsection [O. This opens the possibility 
to investigate the longitudinal expansion of the source |ji 

4.1 Remarks on the fit procedure 

Recently it has been found in the experimental analysis that fits with the YKP para- 
metrisation tend to be more unstable than the corresponding Cartesian fit. We do not 
have the solution of this problem, but we want to present two ideas which might help 
to solve it. 

The YKP parametrisation is constructed in a little bit sophisticated and compli- 
cated way. This provides a relatively simple interpretation of the fit parameters. On the 
other side, the fitting procedure could be more transparent by the usage of a Gaussian 
parametrisation of the following form (and cf. JlT 



C(g,K) - 1 = exp[-fli(K) g i - Ri(K)gf - i#(K)(g°) 2 - 2 R 2 zt {K) m »} . (31) 

This parametrisation might provide better insight to the sources of possible statistical 
uncertainties, since its parameters have clear geometric information concerning the 
Gaussian shape of the correlation function in the (q°, g;)-plane. The YKP parameters 
could be then calculated from the parameters ( |3lD by using the following formulae (for 
Rzt + 0): 



Rl + Rl-^{Rl + RlY-ARl t 
v = , (32) 



Rl-R$ + ^{Rl + Rl?-±R 2 z t Rl - v 2 R 2 t 

R \\ = 2 = 1 + ' (33) 

Rl-Rl + ^Rl + Rtf-ARl t g _ V 2 R 2 
R ° = 2 = 1 + ^ ' (34) 

while R z t = => v = 0, Ru = R z and Rq = Rt- 

The second idea which could lead to the improvement of the fitting is the choice of 
the appropriate reference frame. Although the radii are longitudinally boost invariant, 
the shape of the correlation function is drastically changed when boosting to another 
reference frame. This is shown on Fig. [|. in the YK frame the projection of the 
correlation function into the qi-qo plane is not rotated and has "reasonable" widths. 
When increasing the value of v for fixed R±, Rn, Rq, i.e., when boosting to another 
frame, the projection rotates to the diagonal and becomes very narrow. A shape of 
this kind may cause additional problems in the fitting procedure unless the binning of 
the data is adjusted to the shape of the correlation function. 
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Figure 1: Contour plots of the YKP parametrisation of the correlation 
function projected into the (qi-q°) plane. The values of the radii are: 
i?H = 5 fm, Rq = 2.2 fm/c. The values of the YK velocity on different 
plots are: a) v = 0, b) v = 0.5, c) v = 0.75, d) v = 0.9. 



4.2 Model studies 



To illustrate the main properties of the YKP radii we present here results obtained 
from a model study. The radii are calculated numerically using the formulae dT7|-p!s|) 
and (pC|-23). As a model for the source the emission function of ref.|}4| which is a 
special case of the emission function from ref.p4| is taken. It is given by the formula 



S(x,K) 



M± cosh(?7 — Y) 



exp 



(27t)V27t(At) 



K ■ u(x) 
f 



exp 



(V ~ Vo) 



to 



2R 2 2(Aj?) 2 2(At) 



(35) 

The first part of this emission function describes the geometry of the freeze-out hyper- 
surface. The second is a Lorentz invariant Boltzmann distribution which reflects the 
assumption of local thermal equilibrium at freeze-out. The last exponential describes 
the finite geometrical size of the source in the transverse direction, space-time rapidity 
and proper time. The motion of the different fluid elements of the source is expressed 
by the velocity field u(x). Here we assume longitudinal expansion of the Bjorken type 
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Figure 2: The YK rapidity, denned as Y YK = 0.5 ln(l + v)/(l — v). Transverse 
flow rjf = 0.3, solid lines - pions, dashed lines - kaons. a) Dependence of 
Y YK on the rapidity of the pair Y for fixed K±; b) dependence of Y YK on M± 
for fixed Y; c) dependence of the difference between Y YK and the rapidity of 
the point of maximal emissivity Y LSPS on Y for fixed K±; d) same as c) but 
dependence on Mj_ for fixed Y. 



and a linear transverse expansion profile. Then the components of u(x) are given 
by: 

x y 

u(x) = (cosh?] cosh i]t(r), sinh^(r) — , sinh r]Ar)-, sinhr/ sinh^fr)) (36) 

r r 

with 

r]= : L ]ll t -±l (37) 
' 2 t-z 

and 

r lt{r)= r n f - (38) 

The parameter rjf scales the strength of the transverse flow. 

The temperature T in the calculation was set to 140 MeV, the transverse geometric 
radius R is 3 fm, the space-time rapidity width Ar] = 1.2, the average freeze-out proper 
time To = 3 fm/c, and the mean proper emission duration At = 1 fm/c. Calculations 
were done for pions (m n = 139 MeV) and kaons (mx = 494 MeV). 

The behavior of the YK velocity in this model is plotted in Fig. [2] 0. The results 

2 Figures ^ and ^| are taken from Ref . Jl9| . 
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are presented using the rapidities instead of velocities. The Yano-Koonin rapidity is 
connected with the YK velocity via the usual relation 



Fig. ^a shows the dependence of the YK rapidity on the rapidity of the pair. The linear 
increase is a consequence of the assumed longitudinal expansion with a Bjorken flow 
profile. It is also seen that the lines corresponding to the pairs with higher transverse 
mass are closer to the diagonal. They reach the diagonal in the limit M± — > oo. This 
tendency is evident in Fig. |2|b. The YK rapidity is sensitive mainly to the longitudinal 
rapidity of the pair, the dependence on its transverse mass being weaker. The M± 
dependences for pions and kaons almost coincide; the slight breaking of M±_ scaling is 
caused by the non- vanishing transverse flow [IS]. 



The plots c) and d) of Fig. |2| show the difference between the YK rapidity and the 
rapidity of the point of maximal emissivity of the source for given pair rapidity and 
transverse mass. This difference is appreciable only in the region with non-zero pair 
rapidity and small M±. But already in the region of measurements with kaons the 
difference is very small. 

In Fig. H the characteristic behavior of the YKP radius parameters is plotted. 
The two columns correspond to different physical scenarios: the left one is without 
transverse flow while in the right one the value of r/f is quite big (r/f = 0.6). Before 
embarking on an interpretation of the dependences let us recall that the HBT radii 
measure the homogeneity lengths and not the dimensions of the whole fireball. 

It is clearly seen that the strong transverse flow breaks the common M± scaling 
- the lines for pions and kaons do not coincide. Without the transverse flow the 
only parameter in the emission function is the transverse mass and the M± scaling 
is restored. The influence of transverse flow is most evident in the behavior of R±. 
For i]f = there is no dynamics which could influence the transverse radius, hence 
it is given by the geometric transverse radius of the source. In the real experiment 
there could be another dynamical effects leading to the non-trivial M± dependence, 
e.g. temperature gradients. But the only yet known effect, within the class of thermal 
models, leading to the breaking of common Mj_ scaling is the transverse flow. 

The decrease of Ru is a result of the longitudinal dynamics with strong velocity 
gradient. The fluid element with given z-coordinate moves with the longitudinal ve- 
locity 

v = - f . (40) 

The particle with some longitudinal velocity is effectively produced by the fluid element 
with the same velocity and its surroundings with velocities within the thermal smearing. 
The thermal smearing in this way delimits for higher M± a smaller velocity interval 
than for lower M±. According to eq.(|40|) the longitudinal radius decreases. Again, 
transverse flow destroys the common M± scaling. 

The behavior of Rq is related to the M± dependence of Rn. It is given by the 
geometry of the freeze-out hypersurface. This source freezes out along the hyperbola 
in z-t diagram. Thus events at the freeze-out hypersurface with different longitudinal 
coordinates occur at different times and the longitudinal homogeneity length deter- 
mines partially the effective emission duration (lifetime) ju], ^4|. If the longitudinal 
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Figure 3: Dependences of the YKP radius parameters on Mj_. 
Solid lines - pions, dashed lines - kaons. Left column: rjf = 0, right 

column: rjf = 0.6. The lines for Rq in the case with the transverse 

i 

flow are also compared to the real effective lifetime (t 2 ) 2 . Without 
i 

transverse flow (t 2 ) 2 coincides with Rq. 



homogeneity length vanishes (i.e. at large M±), the lifetime is given by the geometric 
temporal term with width At. 

In the case with transverse flow, the difference between Rq and the effective emission 
duration is plotted too. For high M±, Rq is about 50% bigger than the real lifetime 
of the source, but note, that the assumed value of rjf is rather high. For realistic cases 
we expect rjf about 0.2 - 0.3. If the transverse flow vanishes, Rq measures the lifetime 
exactly. 

5 Conclusions 

The Yano-Koonin-Podgoretskh parametrisation of the two-particle correlation function 
has been presented. Comparing to the older and widely used Cartesian parametrisation 
it provides better insight into the dynamics of the source. 
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A great advantage is the direct measurement of the emission duration provided by 
Ro . The lifetime is the leading term in the model- independent expression for Rq and it 
is measured with some small deviations arising from the non-vanishing transverse flow. 

The fourth fit parameter, the YK velocity or YK rapidity, offers a possibility to 
study the longitudinal dynamics of the source. Its linear increase with the rapidity of 
the pair is a signal for a longitudinal expansion. On the other side, if Y YK does not 
depend on the pair rapidity, this is a signal for a source without longitudinal expansion, 
producing particles with all rapidities from a source at rest in the CM. Please note, 
that these statements are not based on model independent results but on the results 
of model studies within a class of thermal models. 
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